Quantized pumped charge due to surface acoustic waves in a ID channel 
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The adiabatic pumped current through an unbiased one dimensional (ID) channel, connected to 
two ID leads and subject to surface acoustic waves (SAW), is calculated exactly for non-interacting 
electrons. For a broad range of the parameters, quantum interference generates a staircase structure 
of the time-averaged current, similar to experimental observations. This corresponds to integer 
values (in units of electronic charge) of the charge pumped during each period of the SAW. We also 
find staircases for higher harmonics. Quantum interference can thus replace Coulomb blockade in 
explaining the pumped charge quantization, particularly in the SAW experiments. 
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Mesoscopic devices which exhibit quantized dc cur- 
rents may serve as accurate current standards, to be 
used in single-electron metrology. jl| Theoretically, such 
quantization was first suggested by Thouless, j^j for non- 
interacting electrons subject to a slowly moving periodic 
potential which is superimposed on an infinite periodic 
system. This model is an example of adiabatic quan- 
tum charge pumping, ^ a phenomenon which attracted 
much recent theoretical interest. J|-[l8| In such pumping, 
an oscillating potential, with period r = 2tt/u> = 1//, 
generates a time-averaged current I without any bias be- 
tween the two terminals. Under ideal conditions, I has 
the quantized values Afef, with an integer J\f, where e 
is the electron charge. This implies that the charge Q 
transferred during time r is Q = Me. The direction of 
the current is determined by phase shifts between the 
oscillating potential at different locations. 

Quantized pumped currents have been observed in two 
types of experiments. In the first, the periodically mov- 
ing potential is provided by the piezoelectric potential of 
surface acoustic waves (SAW) generated in a quasi-one- 
dimensional (ID) GaAs-AlGaAs channel, |]r^ , po| and the 
observed average acoustoelectric current exhibits steps 
between plateaus at Afef as function of cither the gate 
voltage, V, or the amplitude of the SAW, P. In the sec- 
ond, out-of-phase oscillatory voltages were applied to the 
barriers connecting a quantum dot (QD) to the leads, in 
turnstile-type devices. pT|-p3t In both cases, the appar- 
ent quantization of Q was attributed to the Coulomb 
blockade which dictates the integral number of electrons 
'carried' by the moving potential well. Although such 
explanations may apply to almost closed dots, electron- 
electron interactions should be less important in almost 
open dots (as used e. g. in pl|). 

Adiabatic unbiased quantum pumping for non- 
interacting electrons was recently studied for the 
turnstile-like geometry. [ p"3| When the Fermi energy in 
the leads Ep aligns with the energy of the quasi-bound 
state in the QD, Q/e was found to be close to one. Along 



similar lines, [B the existence of resonant states in a QD 
was found to greatly enhance the magnitude of Q, caus- 
ing it to change sharply as function of some parameters. 
Here we present a simple model for the SAW experiment, 
which exhibits steps and integer plateaus in Q, as func- 
tion of parameters, for non-interacting electrons. Each 
step in Q is accompanied by a narrow peak in the other- 
wise very small time-averaged transmission T. The cal- 
culated staircase plots (e. g. Fig. |l] below), and the ac- 
companying peaks in T, are qualitatively reminiscent of 
the quantum Hall transverse and longitudinal resistances. 
Indeed, both effects follow from steps in Berry-like phases 
of the quantum wave functions, resulting only from quan- 
tum interference. |^,^,|| Interestingly, our results are also 
similar to those observed in the SAW experiments. 

Unlike the turnstile-like case, the piezoelectric poten- 
tial generated by the SAW oscillates with time every- 
where inside the nanostructure, 7isAw( r ,*) = Pcos(u)t — 
q • r) , with the SAW wavevector q, while being heavily 
screened by the 2DEG forming the leads. SAW usually 
satisfy the adiabatic conditions, as hu> is small compared 
to the relevant electronic energy scales. Q In the ab- 
sence of bias, SAW generate a non-zero average current, 
in the direction of q. A realistic treatment of the exper- 
imental geometry jl^] only allowed a calculation at low 
P, yielding Q cx P 2 . The decay of Hsaw m the wide 
banks of the channel is also difficult to treat exactly. In 
order to obtain a simple solvable model, we simplify all 
of these by a ID tight-binding model. The 'channel' is 
made of sites n = 1, 2, TV. It is connected to two 
ID 'leads' via sites 1 and N. The leads, which connect 
to electron reservoirs (with the same chemical potential) , 

are described by Ho = —JJ2 n (\ n )( n + + f° r rl < 
and n > N, with eigenfunctions like e ±lkan and energies 
E = — 2 J cos ka, where a is the lattice constant. The 
'channel' is modeled by 

n osc = ^{e n (*)|rc)(n| - J n (\n)(n + 1| + he)}. (1) 

n 

with J n = Jjj inside the channel, i. e. for 1 < n < N — 1, 
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and Jo = J i, Jn = Jr for the 'contacts' with the leads. 
To model WsawO", i), we choose the site energies inside 
the channel (1 < n < N) as 

e n (t) =V + Pcos[ojt — qa(n — no)], (2) 

where V represents the gate voltage and P > 0, so that 
e n has a maximum (or minimum) in the center of the 
channel n = (N + l)/2 at t = (or r/2). This potential 
acts only inside the channel, imitating the screening of 
the SAW outside. 

Figure [I] shows Q (in units of e) versus V for N = 6 and 
some 'optimal' parameters, at zero temperature (T = 0). 
One clearly observes three plateaus at each end, with 
Q/e very close to ±Af and Af — 0, 1, 2, 3. Generally, 
we observe sharp plateaus up to Af — N/2, with possibly 
several additional rounded peaks or spikes. The steps, at 
Vm, between these plateaus appear to be equidistant; for 
large N we show below that 

W « E F ± (P + 2J D - A(Af + i)), 
A = ga^2PJ D , (3) 

i. e. the steps move outwards (left and right) and 
broaden with increasing P and Jd , as also happens qual- 
itatively in the experiments. (l!]j2(| It would be nice 
to examine this detailed quantitative prediction exper- 
imentally. As one expects, Q — when any one of 
P, Jl, Jd, Qcl, ka or Ep — V vanishes. Indeed, the steps 
become rounded and decrease gradually as P, Ep — V or 
qa approach zero. The rounding begins at the larger A^'s; 
the plateaus at Af = ±1 disappear last. The results re- 
main robust for a wide range of ka, Jl and Jd, provided 
< Jl/J < Jd < P, \E F - V\; Fig. | shows the effects 
of increasing Jp>, and of going to the limit of a completely 
open channel (Jl — 1 and ka = tt/2). In the latter case, 
electron-electron interactions should be unimportant! As 
Eq. (|3j) implies, we also observe steps and plateaus for Q 
versus P at fixed V: at low P, one starts with Q cx P 2 , 
but Q remains very small up to Pq — Ep — V— 2 Jd+A/2. 
Above Pq one observes N/2 steps, at intervals A (which 
now increases with P), and then a gradual decrease to- 
wards zero. Thus, both V and P can be used as triggers 
for on-off switching of the pumped current. 

The 'optimal' graph in Fig. ^ was derived at qa — 
ir/[2(N — 1)], corresponding to a total phase shift ofw/2 
between the end points (n = 1 and N ) , or to a SAW wave 
length A equal to 4L, where L = (N— l)a is the length of 
the channel. Figure || shows results for X/L = 1, 1.5, 2 
and 8; clearly, the staircase deteriorates as A moves away 
from 4L, in both directions. Even as the higher plateaus 
deteriorate, the plateaus at Af = ±1 remain quite robust, 
down to A ~ L as apparently used in the experiments. 
p9[ Note that in the latter, the SAW decays towards the 
ends of the channel, and therefore the relevant L may be 



smaller than estimated, bringing the experiments closer 
to our 'optimal' range. 

Our calculations involve three main steps. In the first, 
we use the adiabatic approximation (neglecting high or- 
ders in time-derivatives of the wave functions) to cal- 
culate the unbiased current. Such calculations usually 
employ the Brower formula which involves integrals 
over the multidimensional space of the time-dependent 
potentials. Here we use an equivalent formula, which is 
more convenient for a general H OS c- If \Xa) denotes the 
instantaneous scattering solution at time t, which results 
from an incoming wave \w~) from lead a (= £, r) with en- 
ergy E, then the instantaneous current from left to right 
is |§ 

where f(E) is the Fermi distribution (which is the same 
in both leads, in the unbiased limit), and the dot denotes 
the time derivative. Below we use T = 0, i. e. E = Ep. 
Using the charge pumped in one period becomes 
Qe = j_^, 2 dtl\. For a symmetric channel (Ji = J r ) 
one has Qi = —Q r = (Qi - Q r )/2 = Q = It. 

The second step involves finding \Xa)- This is similar 
to a static scattering solution: we write 

(n|xl) - A je ikan + A £ e- ikan , n < 0, 

(n\ X \) = B e e ikan , n>N + l, (5) 

with Aqj = 1/ V2Jsin ka (for a unit incoming flux). The 
Schrodinger equations at the sites n = and n = N + 1 
now yield A 0/ e lka + A ee - lka = fa(l)Jg/J and B e e lkNa = 
4>i(N)J r / J ', where {4>e(n)} are the amplitudes of the scat- 
tering solution inside the nanostructure, which obey 

E n >M nn Kt>t(n') = 25 nA i sin kae tka J e A 0/ , 
M n , n , = {g-\E)) ninl =E6 n , n ,-{rl osc ) nnl 

+Sn,n>e ika (5 n ,xJt + 5 n , N J 2 r )/J, (6) 

for 1 < n, n' < N. The solution of these equations is 

4>i{n) = e lka A ^2ismkaJeg n ,i, 

A t = e 2lka A a j[2ismkag 1A Jf/J- 1], 

Bi = eJ- ka[1 - N) A QA 2ismkag NA J g J r /J, (7) 

with similar expressions for cf> r (n), A r and B r . The last 
equation in ([?]) yields the instantaneous ('normal') trans- 
mission, T* = 4| <7jv,i 1 2 ( Ji J r /J) 2 sin 2 ka. In the cases 
of interest here, T is usually very small. However, T* 
has local peaks (as function of cosut) at the N poles of 
g _/v,i 7 i- e. the zeroes of D(coswi) = detM. The time- 
average T = J^J^^dtT* /t thus exhibits peaks wherever 
such poles occur within the period r. 

If TLosc depends on time only via e n , then Eq. (0) yields 



2 



eJ?8mfco r /2 dt J2e n \ g nA \ 2 . (8) 

and shows singularities whenever cos ut comes close 
to one or more zeroes of D within the integration. For 
small (J 2 +J 2 )/ J, these poles have small imaginary parts, 
and Q exhibits large changes whenever cos cot passes near 
such a pole. These steps occur exactly when T has spikes, 
originating from the same poles. 

In the third step, we apply the above adiabatic equa- 
tions to our specific ID channel model. The N x N 
matrix AA = g~ x (E) is now tridiagonal, with A4 n , n — 
E - e„ + e {5nsJi + ^mJ?)/J and M n , n ±i = Jd- 
Since e„ = — coPsm[cot — qa(n — no)], If is equal to co 
times a function of cot, Jaw and Q is independent of 
co. For our symmetric channel (Ji — J r ), D is an order- 
N polynomial in cos cot; we calculate Q by rewriting the 
integrand as a sum over its 2N complex poles in cos cot, 
and using analytic expressions for J Q duu m / '(cos u — z). 

It is interesting to follow these poles as function of the 
various parameters. For this purpose, we show in Fig. 
|] the partial charge Qe(t), resulting from integration of 
Eq. (U) only up to t < t/2, at different values of V. As 
V increases through V\, Qe(t) suddenly exhibits a step 
from zero to one, which appears at t — 0. This step corre- 
sponds to a pole in ij, which enters near cos cot =1. As V 
increases, this step moves to the left. At V = V2, another 
step (from 1 to 2) enters at t = 0. For large N, the time 
interval between two consecutive steps is roughly equal 
to S = qa/oj. As V increases further, both steps move to 
the left, and at V — V3 a step from 2 to 3 enters at t = 0. 
After a narrow intermediate complex state, there begin 
to enter steps of —1, until at V — Ep we observe exactly 
N/2 steps of +1 followed by N/2 steps of —1, ending up 
with Q — 0. A similar build-up of (negative) steps occurs 
when one starts at large positive V, and follows V down 
through V-N, except for the fact that now the new steps 
show up at t — t/2, i. e. coscot = — 1. The physical in- 
terpretation of these results is clear: unlike the Coulomb 
blockade picture, where Af electrons move together from 
left to right, carried by a single minimum of the moving 
potential, in the SAW case treated here Qi changes by 
steps of 1, implying separate motion of individual elec- 
trons, building up to Af after a full period. The picture 
is particularly interesting near V — Ep: during a period, 
Qt(t) exhibits several positive steps, and then an equal 
number of negative steps, ending with no net pumped 
charge. 

The step- like time dependence of Qe{t) immediately 
implies the appearance of higher harmonics (in u>) of 
the pumped current: for V in the A/"'th plateau, we 
approximate 1\ w e'^^L 1 6(t — tj), with tj = t (V) + 
(j — 1)5. The pumped current then becomes 1\ — 
ef{N+ 2 £^=i cos[mw(t -t -(M- l)S/2)]l(m,Af)}, 
with 



l(m,Af) w sin(A/mw<V2)/sin(mw<V2). (9) 

Thus, this m'th harmonic amplitude also exhibits a stair- 
case structure. Although \T(m, 1)| 2 = 1, the plateaus at 
higher Af oscillate with m. It is interesting to note that 
harmonics with m^O survive in the intermediate range 
V ~ Ep: as seen from the center Fig. |], one then has a 
difference like cos[moj(t — ti)] — cos[muj(t — tz)]. It would 
be very interesting to study X(m,Af) experimentally. 

Since the steps in Q(V) arise when poles appear (or 
disappear) when cos tot = ±1, we confirmed that the 
VV's are the solutions of D(V, cos cot = ±1) = 0. Taking 
t = t/2, this equation is equivalent to the set of equations 

(E-V + P cos[(n - n )qa\) 4>(n) = 
-J D [<t>(n+l) + <f>(n-l)] (10) 

for 1 < n < N, with modified equations for n = 1 and 
n = N. Except for the boundaries, these are Mathieu's 
equations, with the symmetric potential having a mini- 
mum at n = no- This can be written as 

-J D [<t>(n + 1) + <f>(n - 1) - 2<f>(n)] 

+P(1 - cos[(n - n )qa])<f)(n) = E(j){n), (11) 

with E = E — V + P + 2Jpj- For our 'optimal' qa = 
ir/[2(N — 1)], one has \n — no\qa < 7r/4, so that the re- 
placement 1 — cos[(n — no)qa] w ^(n — no) 2 (qa) 2 forms an 
excellent approximation. The low lying eigenenergies of 
the isolated channel then correspond to localized states 
around no, which are not very affected by the boundaries. 
For large N, we also use [<fi(n + 1) + <p(n — 1) — 2(f>(n)] ~ 
d 2 <fi/dn 2 . Eq. (|Tl| ) then becomes the Schrodinger equa- 
tion for the harmonic oscillator, yielding the eigenvalues 
Ea/- = A(Af + 1/2), with integer Af, i. e. Eq. (|) with 
the upper sign. The solution for the other sign, associ- 
ated with poles entering at t = 0, follows the replacement 
4>(n) — > (— l) n <j){n). These approximate values agree very 
well with our direct solutions of D(V, ±1) = and with 
the steps in Q, even for relatively small N (e. g. Fig. |l|). 

A few more comments are in place: (a) Although our 
calculated plateaus look like constant integers, in practice 
the function Q(V) remains slightly below Af, reaching a 
smooth maximum around the middle of the 'plateau'. For 
N = 6, the difference (Af - Q(V)) is of order .0001. This 
difference becomes smaller for larger N. (b) Our robust 
plateaus at integer values occur only for the model pre- 
sented here, where the amplitude P has exactly the same 
value for all n in the channel. Modulation of P in space, e. 
g. due to a screening decay, due to a reflected SAW, pc| ] 
or due to random energies {V^}, gives plateaus at non- 
integer values or round the steps. This, and a treatment 
of more complex nanostructures, will be reported else- 
where. p4| ] (c) Our formalism allows for T > 0, where we 
only expect some rounding of the steps, (d) Eq. (|3|) for A 
remains true in the limit N — > 00, a — > 0, (N — l)a = L, 
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when Jdo, 2 — > Ti 2 /(2m*). (e) Finally, we emphasize again 
that all of our results are valid only in the adiabatic limit. 
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FIG. 1. Pumped charge Q (in units of e) versus the 
gate voltage V (in units of J) for N = 6, P = 8J, 
Jd — J, Ji = Jr = Jl — -4J, qa = 7r/10, ka = 7r/100. 
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FIG. 2. Same as Fig. |, but with J D = 3J (left) and 
Jl = J, ka = 7r/2 (right). 
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FIG. 3. Same as Fig. |, but with qa = 2tt/5, 4tt/15, tt/5 
and 7r/20, corresponding to \/L = 1, 1.5, 2 and 8. 
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FIG. 4. Partial pumped charges Qe, in units of e, up to 
time t within a period, for the same parameters as in Fig. [l|, 
with V/J — —8.6, — 2 and 5.3. The i-axis shows cut between 
— 7r and 7T. 
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